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Abstract 

Using the leading twist approach to nuclear shadowing, which is based on the relationship be- 
tween nuclear shadowing and diffraction on a nucleon, we calculate next-to-leading order nuclear 
parton distribution functions (nPDFs) and structure functions in the region 0.2 > x > 1CT 5 and 
Q 2 > 4 GeV 2 . The uncertainties of our predictions due the uncertainties of the experimental input 
and the theory are quantified. We determine the relative role of the small (~ Q 2 ) and large (3> Q 2 ) 
diffractive masses in nuclear shadowing as a function of x and find that the large mass contribu- 
tion, which is an analog of the triple Pomeron exchange, becomes significant only for x < 10 . 
Comparing our predictions to the available fixed-target nuclear DIS data, we argue, based on the 
current experimental studies of the leading twist diffraction, that the data at moderately small 
x ~ 0.01 and Q 2 ~ 2 GeV 2 could contain significant higher twist effects hindering the extraction 
of nPDFs from that data. Also, we find that the next-to-leading order effects in nuclear shad- 
owing in the ratio of the nucleus to nucleon structure functions F2 are quite sizable. Within the 
same formalism, we also present results for the impact parameter dependence of nPDFs. We also 
address the problem of extracting of the neutron F2 n (x,Q 2 ) from the deuteron and proton data. 
We suggest a simple and nearly model-independent procedure of correcting for nuclear shadowing 
effects using F 2 A /F 2 D ratios. 

PACS numbers: 24.85.+p,13.60.Hb 
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I. INTRODUCTION 



One way to analyze the microscopic structure of atomic nuclei is to study the distribu- 
tion of quarks and gluons, as well as their correlations, in nuclei. These nuclear parton 
distribution functions (nPDFs) can be accessed using various deep inelastic scattering (DIS) 
processes: Inclusive scattering of leptons, high-mass dimuon production using proton beams, 
exclusive electroproduction of vector mesons. None of the above processes determines nPDFs 
comprehensively, only taken together do these experiments provide stringent constraints on 
nPDFs. 

The discussion of the present paper is centered around the nuclear effects of nuclear 
shadowing and antishadowing (enhancement), which affect nPDFs at small values of Bjorken 
variable x, 10~ 5 < x < 0.2. Nuclear shadowing of nPDFs is developing into an increasingly 
important subject because it is involved in the interpretation of the RHIC data on jet 
production, evaluation of hard phenomena in proton-nucleus and nucleus-nucleus collisions 
at the LHC, estimates of the black limit scattering regime in DIS, etc. 

The major obstacle that hinders our deeper knowledge of nPDFs at small x is that, 
up to the present day, all experiments aiming to study nPDFs are performed with fixed 
(stationary) nuclear targets. In these data, the values of x and Q 2 are strongly correlated 
and one measures nPDFs essentially along a curve in the x — Q 2 plane rather than exploring 
the entire plane. Moreover, for Q 2 > 1 GeV 2 , the data cover the region x > 5 x 10 -3 , where 
the effect of nuclear shadowing is just setting in. As a result, when one attempts to globally 
fit the available data by modeling nPDFs at some initial scale Qq and then performing QCD 
evolution, various groups jl} Q, S U, 0] produce significantly different results. 

An alternative to the fitting to the data is to combine the Gribov theory p, which 
relates the nuclear shadowing correction to the total hadron-deuteron cross section to the 
cross section of diffraction off a free nucleon, with the Collins factorization theorem ?J 
for hard diffraction in DIS. The resulting leading twist theory of nuclear shadowing was 
developed in $\ and later elaborated on in ^]. 

The Gribov theory has been applied to the description of nuclear shadowing for many 
years. First it was done in the region of small Q 2 , where generalized vector dominance 
model gives a good description of diffraction, see review in 10J, and later in deep inelastic 
region, where large diffractive masses donate □ A o^ofsuccessM^e. 
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calculations were performed 



12 



13 



14 



15 



16 



17 



18, 



19. 



before the experimental data 



from HERA became available. A calculation constrained to reproduce the HERA data using 
the Gribov theory was presented in j^]. It focuses on the calculation of nuclear shadowing 
for F^ at intermediate Q 2 where leading and higher twist effects are equally important. A 
fair agreement of the data with the Gribov theory has been found. However, this approach 
does not involve the use of the Collins factorization theorem and, hence, does not address 
nPDFs (see a detailed comparison in Sect. II V|) . 

The present work extends the calculation of nPDFs of [l| with an emphasis on the the- 
oretical ambiguity and accuracy of the predictions and makes a comparison to fixed-target 
nuclear DIS data. In particular, we demonstrate that 

• The theory of leading twist nuclear shadowing and QCD analysis of hard diffraction 
at HERA enable one to predict in a model-independent way the next-to-leading order 
nPDFs for 10~ 5 < x < 10~ 2 with 30% accuracy, Fig. El For larger x, 10" 2 < x < 
0.1 — 0.2, there appears an additional effect of nuclear antishadowing that requires 
modeling and whose uncertainty is larger. In addition, the HERA diffractive data 
for xjp > 0.01 contains a sub-leading Reggeon contribution, which adds additional 
ambiguity to our predictions, especially for x > 0.01. 

• The interactions with N > 3 nucleons (which is a model-dependent element of the 
Gribov approach) give negligible contribution in the NMC fixed-target nuclear DIS 
kinematics, see Fig. |HJ The A-dependence of the NMC data for x ~ 0.01 is reasonably 
well reproduced, see Fig. H31 

• The failure to describe the absolute value of the F^ / (AF^) ratios of the available fixed- 
target data for 0.003 < x < 0.02 and Q 2 < 3 GeV 2 likely indicates the presence of 
significant higher twist effects in the data. Indeed, when the leading twist shadowing 
is complemented by higher twist effects, which are modeled by p, <fi and to meson 
contributions in the spirit of vector meson dominance, the agreement with the data 
becomes fairly good, see Figs. El an d ED All this signals that any leading twist 
QCD analysis of the available data is unreliable for 0.003 < x < 0.02. 

• The next-to-leading order (NLO) effects in the F% /(AF^) ratios are found to be quite 
sizable. This means that it is not self-consistent to use the leading order parameteri- 
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FIG. 1: Gribov's theorem The forward hadron-deuteron rescattering amplitude, which gives 
rise to nuclear shadowing, is proportional to the differential hadron-nucleon diffractive cross section 
at t ~ 0. 

zations of nPDFs in the NLO QCD calculations, see Fig. ED 

In short, the main goals of the paper are to give a concise summary of the leading 
twist theory of nuclear shadowing, to assess the theoretical uncertainties of the resulting 
predictions and to make a comparison to nuclear DIS data. We attempted to give a self- 
contained presentation and, hence, this paper can be rightfully considered as a guide to 
leading twist nuclear shadowing. 



II. LEADING TWIST THEORY OF NUCLEAR SHADOWING 

In this section, we review the leading twist approach to nuclear shadowing developed in 
and further elaborated on in jjj. 

The approach is based on the 1969 work by V. Gribov |6j, where the following theo- 
rem was proven. Let us consider hadron-deuteron scattering at high energies within the 
approximation that the radius of the deuteron is much larger than the range of the strong 
interaction. Then the shadowing correction to the total cross section is expressed in terms 
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of the differential diffractive hadron-nucleon cross section. This is demonstrated in Fig. 
The forward hadron-deuteron rescattering amplitude giving rise to the nuclear shadowing 
correction contains the hadron-nucleon diffractive amplitude (denoted by the shaded blob) 
squared. 

The relationship between nuclear shadowing and diffraction was used in the analysis 
of parton densities of deuterium and other nuclei by Frankfurt and Strikman in For 
deuterium and other sufficiently light (low nuclear density) nuclei, nuclear shadowing and 
diffraction on the nucleon are related in a model-independent way using the Gribov theo- 
rem j^2|. 

The generalization to heavy nuclei involves certain modeling of multiple rescattering con- 



tributions, which, however, is under control 



23j | . Below we shall recapitulate the derivation 



of the leading twist nuclear shadowing for nPDFs, which can be carried out in three steps. 
Step 1. The shadowing correction arising from the coherent interaction with any two nu- 

(2) 

cleons of the nuclear target with the atomic mass number A, SF^ (the superscript (2) serves 
as a reminder that only the interaction with two nucleons is accounted for), is expressed in 
terms of the proton diffractive structure function F 2 D ^ (the superscript (4) indicates the 
dependence on four kinematic variables) as a result of the generalization of the Gribov result 
for deuterium (see also Ref. This does not require decomposition over twists and is 

(2) 

therefore valid even for the case of real photon interactions. The shadowing correction SF^ A 



reads 



6F%(x,Q 2 ) = A(A ~ 1] lGnTle 



i-v) 2 f f°° f°° f XF '° 

d 2 b / dz\ l dz2 / dxjp 



vj 2 



Z\ 



X F 2 D{4) (/3, Q 2 , Xjp, t) \ t=tininPA (b, Z l ) PA (b, z^rn^-z,) 



(1) 



with r) the ratio of the real to imaginary parts of the diffractive scattering amplitude; z±, 
Z2 and b the longitudinal (in the direction of the incoming virtual photon) and transverse 
coordinates of the nucleons involved (defined with respect to the nuclear center); (3, xp 
and t the usual kinematic variables used in diffraction. Throughout this work, we use 
(3 = x/xp. Equation (JTJ) uses the fact that the t-dependence of the elementary diffractive 
amplitude is much weaker than that of the nuclear wave function, and, hence, F^(4) can 
be approximately evaluated at t — t min ~ 0. All information about the nucleus is encoded 
in the nucleon distributions pA(b, Zi), see Appendix 1X1 for details. Finally, xjp$ is a cut-off 
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parameter (xp t o = 0.1 for quarks and xp^ = 0.03 for gluons), which will be discussed later 
in the text. 

The origin of all factors in Eq. (Q) can be readily seen by considering the corresponding 
forward double rescattering Feynman diagram (see Fig.|2I), which accounts for the diffractive 
production of intermediate hadronic states by the incoming virtual photon: 

• The combinatoric factor A(A — l)/2 is the number of the pairs of nucleons involved 
in the rescattering process. 

• The factor 16tt provides the correct translation of the differential diffractive to the 
total rescattering cross section (see the definition later), as required by the Glauber 



theory 
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The factor (1 — irj) 2 /(l + rj 2 ) is a correction for the real part of the diffractive scattering 
amplitude A. Since the shadowing correction is proportional to (ImA) 2 , while the total 
diffractive cross section is proportional to \A\ 2 , the factor (1 — it]) 2 /{I + rf ) emerges 
naturally, when one expresses nuclear shadowing in terms of the total diffractive cross 
section (diffractive structure function). 

The integration over the positions of the nucleons is the same as in the Glauber theory. 
Similarly, because the recoil of the nucleons is neglected (the transverse radius of the 
elementary strong amplitude is much smaller than the scale of the variation of the 
nuclear density), both involved nucleons have the same transverse coordinate b. 

The integration over xjp represents the sum over the masses of the diffractively pro- 
duced intermediate states. 

In order to contribute to nuclear shadowing (not to break up the nucleus in its tran- 
sition from the |in)-state to the (out|-state), the virtual photon should interact with 
the nucleons diffractively. The product of the two diffractive amplitudes (depicted as 
shaded blobs in Fig. |2J) gives the diffractive structure function of the proton F^ 4 \ 
Also note that we do not distinguish between diffraction on the proton and neutron 
in the present work, as the corresponding diffractive amplitudes are equal at small x. 

The effect of the nucleus is given by the nucleon densities pA{b, zA. For the sufficiently 
heavy nuclei that we consider, nucleon-nucleon correlations can be neglected and the 




FIG. 2: The forward 7*-nucleus rescattering amplitude that gives the principal contribution to 
nuclear shadowing. 

nuclear wave function squared can be approximated well by the product of individual 
PA(b,Zi) for each nucleon (the so-called independent particle approximation). 

• The factor e lx P m N{zi-z 2 ) j s a consequence of the propagation of the diffractively pro- 
duced intermediate state between the two nucleons involved. 



H and hard diffractive DIS Q 



Step 2. The QCD factorization theorems for inclusive |25( and hard diffractive DIS 
can be used to relate the structure functions in Eq. (JTJ to the corresponding - inclusive and 
diffractive - parton distribution functions. Since the coefficient functions (hard scattering 
parts) are the same for both inclusive and diffractive structure functions, the relation between 
the shadowing correction to nPDFs and the proton diffractive parton distribution functions 
(PDFs) is given by an equation similar to Eq. (0). The shadowing correction to the nPDF 
of flavor j, fj/A-, fifj/A' * s elated to the proton (nucleon) diffractive PDF of the same 
flavor 



^(x,g 2 ) = ^^i6^ e 



1 + v 



' d 2 b l dz\ \ dz 2 / dxjp 



2 



oo J zi 



x/^ 4) (A Q 2 , xjp, t)\ t=tminPA (b, z 1 )p A (b, Z2 ) e ^m N (^ 2 ) 



(2) 
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Equation (j2J) is very essential in several ways. Firstly, it enables one to evaluate nuclear 
shadowing for each parton flavor j separately. Secondly, since the diffractive PDFs obey 

(2) 

leading twist QCD evolution, so does the shadowing correction of^,' A . This explains why 
the considered theory can be legitimately called the leading twist approach. Since Eq. (J2J) 
is based on the QCD factorization theorem, it is valid to all orders in a s . Hence, if ijyjy 
is known with the next-to-leading order (NLO) accuracy, as is the case for the used HI 
parameterization for we can readily make predictions for NLO nPDFs. 

Step 3. Equation (J2J) is derived in the approximation of the low nuclear thickness, and 
it takes into account only the interaction with two nucleons of the target. The effect of 
the rescattering on three and more nucleons can be taken into account by introducing the 
attenuation factor T(b, z±, z 2 ) (see for example [lcl]). 

T{b, z h z 2 ) = e-^X 1 -^ L? dz PA {b, z) ^ (3) 

where the meaning of a 3 eS should become clear after the following discussion. Let us consider 
sufficiently small values of Bjorken variable x such that the factor e lx P m N(zi-z2) j n q 
can be neglected. Then, introducing ai s as 



^•^^otttt)! ^aw,^)u te , (4> 

Eq. (J2J) can be written in the form equivalent to the usual Glauber approximation 
8f-%(x,Q 2 ) « A[A ~ l \ l-r, 2 y cS {x^)f j/N {x,Q 2 ) 



X 



//*OG /'OC 
d 2 b j dz x j dz 2 p A {b,z 1 )p A {b,z 2 ) , (5) 



where fjm is the proton inclusive PDF. Therefore, it is clear that thus introduced a J eS has the 
meaning of the rescattering cross section, which determines the amount of nuclear shadowing 
in the approximation of Eq. (jSJ). Hence, it is natural to assume that the same cross section 
describes rescattering with the interaction with three and more nucleons, as postulated by 
the definition of the attenuation factor T(b, zi,z 2 ) by Eq. ©. In the language of Feynman 
diagrams, the assumed form of the attenuation factor implies that the diffractively produced 
intermediate state rescatters without a significant change of mass with the same cross section 
on all remaining nucleons of the target, as depicted in Fig. El for the case of the triple 
scattering. The approximation of elastic rescattering is not important at small enough x 
where longitudinal distances are much larger than the nuclear size, see discussion below. 




FIG. 3: The forward 7*-nucleus triple scattering amplitude. 

Corrections to the elastic rescattering approximation can be estimated by taking into 
account the effects of fluctuations of the strength of the rescattering interaction. Modeling 
of these effects was performed in [23] with the conclusion that for a wide range of cross 
section fluctuations, the reduction of nuclear shadowing (for fixed cr e g ) remains a rather 
small correction for all nuclei. 

After introducing the attenuation factor into Eq. (J2J), the complete expression for the 
shadowing correction, dfj/A, becomes 



[l-ir]) 2 f , 2 /' x mru 



6f j/A (x,Q 2 ) = A{A ^ lforfte , 
2 [ 1 + T) z 

x/^ 4) (A Q 2 , x P , t min ) PA (b, Zl ) PA (b, 22 ) e ™(---) e - (A/2)(1 -^« dzpA ^ z) 



d b I dzi I dz 2 / dxjp 

— OO J Zl Jx 



I (6) 

This is our master equation (see also Eq. (fTf|0 . It contains several sources of model- 
dependence and theoretical ambiguity. First, the attenuation factor T(6, z\, z-i) assumes 



that multiple rescatterings can be described by a single rescattering cross section 59] a 3 eS , 
i.e. cross section fluctuations are neglected in the interaction with three and more nucleons. 
Note that in the phenomenologically important kinematic region of fixed-target experiments, 
x > 0.01 and Q 2 > 2 GeV 2 , the uncertainty associated with the attenuation factor T(b, Z\, z^) 
is negligible since the rescattering contribution to shadowing is small, see Fig. |H| Second, the 
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ution scale we take Ql = 4 GeV 2 : This is the lowest value of Q 2 of the HI diffractive 
2t|. Nuclear PDFs at Q 2 > Ql are obtained using the NLO QCD evolution equations. 



necessity to introduce the parameter xjp^ is a consequence of the fact that Eq. (jHJ) applies 
only to the region of nuclear shadowing: The transition to the region of the enhancement of 
nPDFs should be modeled separately. This is the role of the parameter xp$. Third, there 
are experimental uncertainties in the determination of the diffractive PDFs which we 
use as an input in Eq. (0). 

Equation © defines the input nPDFs for the DGLAP evolution equations. As a starting 
evo 
fit 

Therefore, we predict that nuclear shadowing is a leading twist phenomenon. 

In the small- x limit, which for practical purposes means x < 10 -3 , the factor e lx P m N{zi-z 2 ) 
in Eq. © can be safely omitted, which results in a significant simplification of the master 
formula (after integration by parts two times) 

5fM^Q 2 ) = 2{1 ' 1/A l/ MX,Q2) Re (/ d 2 b(e-^_ 1 + LT{b)) ^ , (7) 

where L = A/2 (1 - ir}) a J eS ; T(b) = f^dz p A (b, z). 
In the heavy nucleus limit (A — > oo) and at fixed a e g, 

fj/A(x,Q 2 ) 5f j/A (x,Q 2 ) _ 2nR 2 A 

Af j/N (x,Q 2 ) Af, /N (x,Q 2 ) Aai s ' U 

where Ra is the nuclear size. As can be seen, in the theoretical limit of infinitely heavy 
nucleus, nuclear shadowing equals the ratios of the nuclear to nucleon sizes, i.e. it is a 
purely geometrical effect. At the same time, in the Q 2 = const, W — > oo limit, when the 
leading twist approximation is violated and the radius of the strong interaction becomes 
larger than Ra, the a r *A/&-y*N ratio should approach unity 26^ . 

III. PARAMETERS AND UNCERTAINTIES OF THE METHOD 

The master equation © uses as input the information on hard diffraction in DIS on the 

r, n 

proton, which was measured at HERA by ZEUS |27[ and HI 28] collaborations. We use the 
HI parameterization of 13 (note the superscript (3) indicating that the t-dependence 
of diffraction is not measured), which is based on the QCD analysis of the 1994 HI data 

D: 



|28| (we use Fit B, see Appendix A of Ref. |9j). The choice of the HI parameterization is 
motivated by the following observations: 
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It is available in an easily accessible and usable form, see 
Ref. Q. 



and also Appendix A of 



The diffractive jet production in DIS at HERA data 30] is best described by the HI 
parameterization. The fit of Alvero, Collins, Terron and Whitmore |3l| somewhat 
overestimates the data. Another parameterization available in the literature, that of 



Hautmann, Kunszt and Soper 
diffractive data. 



321 ] . is not based on the detailed fit to the available 



• The 1994 HI fit is in a fair agreement with the most recent 1997 HI data |33j. However, 
the 1997 HI data indicates that the gluon distribution of the 1994 fit is too large by 
about 25%. Hence, in our analysis we multiplied the gluon diffractive distribution of 
H by 0.75. 

Since the diffractive PDF f^ffi enters Eq. (JSj) at t ~ 0, one has to assume a certain t- 
dependence in order to be able to use the HI results for the t-integrated ' . The common 
choice is to assume that 

f^\p,Q 2 ,x ]P ,t)=e B > t f^\f3,Q 2 ,x ]P ,t^0), (9) 

so that after the integration over t, one obtains 

f^\(3,Q 2 ,xjp } t^ 0) = Brfffitf&xr) , (10) 

where Bj is the slope of the t-dependence of jfi/jy ■ A priory there is no reason why the 
slope Bj should be equal for all parton flavors j and, hence, we introduce its explicit flavor 
dependence. In our analysis we use the following values for Bj. For all quark flavors, we 
use B q = 7.2 ± l.l(stat.)^Q;g(syst.) GeV~ 2 , which is determined by the measurement of the 
t-dependence of the diffractive structure function F®^\ as measured by the ZEUS collab- 
oration [34]. Of course, the diffractive slope should increase with decreasing x (diffractive 
cone shrinkage). However, since the experimental error of the value of B q is large and no 
measurements of the ^-dependence of B q are available, any theoretically expected logarith- 
mic increase of B q will be within the quoted experimental errors. Hence, it is sufficient to 
use the x-independent B q . 

The slope of the gluon PDF, B g , could be different from B q . If the gluon-induced diffrac- 
tion is dominated by small-size (compared to typical soft physics sizes) partonic configura- 
tions in the projectile, B g could as low as the slope of J/ip diffractive production measured 

12 



at HERA which is substantially lower than B q . To reflect the uncertainties in the value of 
B g we examined two scenarios: B g = 4 + 0.21n(10~ 3 /a;) GeV -2 and B g = 6 + 0.25 ln(10~ 3 /x) 
GeV -2 . The first one corresponds to the lower end of the values of the J/ip photoproduction 
slope reported at HERA 35j, while the second one is close to B g and to the J/ip slope 
reported in j36j. 

The analysis of recent ZEUS data on the slope of j^l reports the values similar 

to those reported in 3^|. The analysis also provides information on the xjp-dependence of 
the slope though not on the /3-dependence of the slope (average (3 for the data sample is 
growing with a decrease of xjp). Overall, the data appear to be consistent with the Regge 
factorization and for the most of the xjp range, the gluon diffractive PDF appears to give a 
significant, if not the dominant, contribution. This suggests that our model with a higher 
value of B g is closer to the data, though in view of the lack of the data on the slope of the 
gluon induced diffraction at large (5, we feel necessary to keep the lower B g model as well. 

The analysis of the 1994 HI data j3] showed that at large xjp, the successful fit to the 
data requires both the Pomeron and the Reggeon contributions. In our numerical analysis 
we include only the dominant Pomeron part, because the subleading Reggeon contribution 
begins to play a role only for x > 0.01, where the theoretical ambiguities are large anyway. 
Therefore, in our analysis we use the Reggeon contribution only to estimate its contribution 
to the overall uncertainty of our predictions, see Appendix |B] for details. 

Using Eq. (JTHJ), the rescattering cross section a J eS becomes 

where rj is the ratio of the real to imaginary parts of the diffractive amplitude A. This 
ratio can be related to the intercept of the effective Pomeron trajectory, ajp(0), using the 
Gribov-Migdal result Q 

»7«|(a^(0)-l) = 0.32, (12) 

where the HI value for ajp(0) was employed. 

The results of the evaluation of a ] eS are presented in Fig. 0J The left panel presents cr e g 
for anti u-quarks and the right panel is for the gluons, both cases for Q 2 = 4 GeV 2 . The 
error bands around the central curves represent the uncertainty in the determination of a c g. 
This uncertainty comes from the uncertainties in B q , f®^\[3,Q 2 } xip) (taken to be 25%) 
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FIG. 4: The effective cross section a e s for the anti u-quark and gluon channels at Qq = 4 GeV 2 . 
The error bands represent the uncertainty in the predictions discussed in the text. 

and the choice of xpp added in quadrature. Two solid curves for the gluon case correspond 
to the two scenarios for the slope B g discussed above. 

Now we would like to examine which values of the diffractive masses or (3 contribute to a e s ■ 
At very high energies (small x), one enters the regime analogous to the triple Pomeron limit 
of hadronic physics, which corresponds to (3 = Q 2 /{Q 2 + M\) <C 1. In this case, one may 
need to resum logs of energy in the diffractive block (logs of (3). However, deviations from 
DGLAP are expected only for (3 < 10 -3 , which is beyond the x range that we consider |39|- 
At extremely small x, the contribution of small (3 may become dominant. This contribution 
was evaluated within Color Condensate model in |40( neglecting the large (3 contribution. 
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FIG. 5: The ratio R at Qq = 4 GeV 2 . The solid curves correspond to /3 max = 0.5; the dashed 
curves correspond to /? max = 0.1; the dotted curves correspond to /3 max = 0.01; the dot-dashed 
curves correspond to /3 max = 0.001. 

To analyze at what x small (3 become dominant, it is convenient to introduce the ratio 
R defined as follows 



-R(Anax, X ) - rXlDn , „D(?C\,„ — , ■ ( 13 ) 



The ratio R for w-quark and gluon channels at Ql = 4 GeV 2 is presented in Fig. 03 In 
the figure, the solid curves correspond to /3 max = 0.5; the dashed curves correspond to 
An ax = 0.1; the dotted curves correspond to /5 max = 0.01; the dot-dashed curves correspond 
to /3 max = 0.001. 
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From Fig. |5] one can see how much different /3-regions contribute to nuclear shadowing. 
For instance, taking x = 10 -3 , which roughly corresponds to the smallest x which could be 
reached in the RHIC kinematics, one sees that large diffractive masses that correspond to 
(3 < 0.1 (dashed curve) contribute 20% to nuclear shadowing in the quark channel and 30% 
to nuclear shadowing in the gluon channel. Therefore, Fig.Elindicates that if the Color Glass 
Condensate model is implemented in a way consistent with the HERA diffractive data, it 
predicts a very small fraction of total shadowing for the RHIC kinematic range. 

For completeness, we rewrite our master equation, Eq. (JUJ), in the form which explicitly 
includes the diffractive slope Bj 



2 I + 77 

X/^ 3 ^, Q\ X P ) PA (b, Zl ) PA (b, Z2 ) e ™Frn N{zi -z 2 ) e -(A/2)(l-ivK a ft dzp A (b,z) 



poo rxpfi 

d z b / dz\ I dz 2 / dxjp 

CO J Z\ J X 



(14) 



Ij/N 

We would like to point out that while the leading twist theory of nuclear shadowing is 
applicable to the partons of all flavors, see Eq. (j!4j). using the low-x HERA diffractive data, 
which is heavily dominated by the Pomeron contribution, we cannot make any quantitative 
predictions for nuclear shadowing of the valence quarks in nuclei. Nuclear shadowing for 
the valence quarks is driven by the t-channel exchanges with non- vacuum quantum numbers 
(Reggeon contribution), whose contribution is largely lost in the kinematic region of the 
HERA data. In practical terms, this means that Eq. (|14|) should be applied to evaluate 
nuclear shadowing for the antiquarks and gluons only. 

As mentioned above, Eq. (|14|) cannot describe nuclear modifications of PDFs at x > 0.1 
for the quarks and x > 0.03 for the gluons, where nuclear antishadowing and the EMC 
effects dominate. For a comprehensive picture of nuclear modification for all values of x, we 
refer the reader to the review in ^J. However, since we use Eq. (|14jl to evaluate nPDFs at 
some input scale for QCD evolution, in order to provide sensible results after the evolution, 
we should have a reasonable estimate of nPDFs for all x. We adopted the picture of nuclear 
modification of PDFs developed in [ill . which suggests that antiquarks in nuclei are 
not enhanced and the gluons are antishadowed and which uses the constraints based on the 
baryon and energy- momentum conservation sum rules. In our case, like in jjj, we model 
the enhancement of the gluon nPDF in the interval 0.03 < x < 0.2 with a simple function 
a(0.2 — x)(x — 0.03) and choose the free coefficient a by requiring the conservation of the 
momentum sum rule for nPDFs. For instance for 40 Ca, this requirement gives a ~ 30 and 

16 



about 2-3% enhancement of the fraction of the total momentum of the nucleus carried by 
the gluons, in accord with the analysis of 



IV. LEADING TWIST NPDFS AND STRUCTURE FUNCTIONS 

The master equation ()14|) allows one to determine NLO nPDFs at the input scale Qq = 4 
GeV 2 . As an example of such a calculation, we present ratios of the nuclear (Ca-40) to 
free proton PDFs and the ratio of the nuclear to the free nucleon structure function, = 
(F$ + -F^)/2, at Q 2 = 4 GeV 2 by solid curves in Fig. El The shaded error bands around 
the solid curves indicate the uncertainty of the predictions. For comparison, LO predictions 
of Eskola et al. [1] (based on the LO fit to the DIS and Drell-Yan nuclear data) for the 
corresponding ratios are given by the dashed curves. 

Note that for the gluon ratio, we give two predictions corresponding to two versions of 
the diffractive slope B g discussed earlier. Also, since we do not predict nuclear shadowing 
for the valence quarks, this information should be taken from elsewhere. In our analysis, we 
use the parameterization by Eskola et al. |l| (see the upper left panel). 



421. 



For the parameterization of the proton PDFs, we used the NLO fit CTEQ5M 
Figure [3 presents the Q 2 -evolution of the ratios in Fig. El The solid curves correspond to 
Q 2 = 4 GeV 2 ; the dashed curves correspond to Q 2 = 10 GeV 2 ; the dot-dashed curves corre- 
spond to Q 2 = 100 GeV 2 . The leading twist character of the predicted nuclear shadowing 
is apparent from this figure: The shadowing correction decreases slowly with increasing Q 2 
and there is still rather significant nuclear shadowing at Q 2 = 100 GeV 2 . 

One can see from Fig.EJthat our predictions at the lowest values of Bjorken x significantly 
differ from those by Eskola et al. However, one should keep in mind that we make our 
predictions to the NLO accuracy, while the fitting to the nuclear DIS data in |l| is done 
to the LO accuracy. One should also note that the parameterization of Eskola et al. 
assumes that at small x, the ratios / (AF^) and g^/ {Ag^) become equal and stay constant 
(saturate). We point out that: 

• Figure El presents our predictions for the shapes of the nPDFS for 40 Ca, which are to 
be used as an input for QCD evolution at the scale Qo = 2 GeV. This choice of Qo 
is motivated by the fact that the 1994 HI diffractive data has Q 2 > 4.5 GeV 2 and 
the QCD fit to the data of j3| starts at Q 2 = 3 GeV 2 . Results of such evolution are 
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FIG. 6: The ratio of nuclear to proton NLO parton distributions and the nuclear to free nucleon 
inclusive structure functions F2 in 40 Ca at Q = 2 GeV. The leading twist theory results (solid 
curves and the corresponded shaded error bands) are compared to the LO predictions by Eskola 
et al. |jl (dashed curves). 

presented in Fig. 

• Leading twist theory predicts much more significant nuclear shadowing for quarks and 
gluons than the fits to the fixed-target data of Eskola et al. [lj. The latter assumed 
that shadowing saturates for small x, e.g. for x < 3 x 1CT 3 for 40 Ca, and that higher 
twist effects are negligible. 

• Nuclear shadowing for the gluons is larger than for the quarks due to the dominance 
of gluons in diffractive pdfs (this was further confirmed by the recent ZEUS data (37 1. 
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FIG. 7: Scaling violations for nuclear PDFs in 40 Ca. The solid curves correspond to Q 2 = 4 GeV 2 ; 
the dashed curves correspond to Q 2 = 10 GeV 2 ; the dot-dashed curves correspond to Q 2 = 100 
GeV 2 . 

• Within our model of antishadowing for the gluons by the simple function a(0.2 — x)(x— 
0.03), significant variations of the parameter a still lead to the conservation (with 
accuracy better than 1%) of the parton momentum sum rule. Hence, the amount of 
antishadowing for the gluons is not sensitive to the low-x behavior of the gluons. 



Should we compare our predictions to those 
the updated fit by Hirai, Kumano and Nagai 



3 by 



3y Hirai, Kumano and Miyama or to 
the disagreement in the shadowing 



predictions, especially for the gluons, would be much larger. For the comparison of 
the parameterization* of Q and the first of Q one ean eonsult Q. 

Comparing to the parameterization suggested in the work of Li and Wang [sj , we again 
find a strong disagreement in the quark channel and, a surprisingly good agreement 
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for the gluons. However, the parameterization of B are not based on the detailed 
comparison to all available fixed-target data, but rather on the need to fit the RHIC 
data within the HIJING model. Note also that this parameterization does not include 
the enhancement either in the quark or gluon channels and, hence, violates the exact 
QCD momentum sum rule for the parton densities. 

The only NLO QCD fit to nuclear DIS data by de Florian and Sassot [5] produces 
very small nuclear shadowing for nPDFs, which is inconsistent with our predictions as 
well as with the predictions of jllQ, sj]- 

Our predictions for nPDFs and the structure function F£ for the nuclei of 12 C, 40 Ca, 
110 Pd, 197 Au and 206 Pb and for the kinematic range 10~ 5 < x < 1 and 4 < Q 2 < 10, 000 
GeV 2 have been tabulated. They are available in the form of a simple Fortran program 
from V. Guzey upon request, vadim.guzey@tp2.rub.de, or from V. Guzey's web page, 



http:/ /www.tp2. rub.de/~vadimg/index. html 



We also study the importance of the effect of multiple rescatterings, which is described by 
the attenuation factor T(b) in Eq. (fT4*j) . Figure |H1 compares the result of the full calculation 
of the ua/{Aun) ratio for 40 Ca at Q = 2 GeV (solid curve) with the calculation, when the 
rescattering effect was ignored (dashed curve), i.e. T{b) was set to one in Eq. (fl"4j) . As seen 
from Fig. [HI the rescattering effect becomes unimportant for x > 0.005, i.e. in the kinematics 
of the fixed-target nuclear DIS experiments. 

Nuclear shadowing corrections to nPDFs become significantly larger, when one considers 
the interactions with the target nucleus at small impact parameters. Indeed, since the 
density of nucleons in the center of the nucleus is larger than the average nucleon density, 
choosing small impact parameters corresponds to the increase of the number of scattering 
centers. Introducing the impact parameter dependent nPDFs, fj/ A (x, Q 2 ,b), as was done in 

B 

J d 2 bf j/A (x, Q 2 , b) = f j/A (x, Q 2 ) , (15) 

the nuclear shadowing correction to the impact parameter dependent nPDFs can be readily 
found from Eq. ()14j) by simply removing the integration over the impact parameter b 



Sf j/A (x,Q 2 ,b) = A{ \ — 16 BjiiTZe 



2 roo poo rxp.o 



dzi J dz 2 I '/■'•//• 
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FIG. 8: The ua/(Aun) ratio for 40 Ca at Q = 2 GeV. The solid curve is the result of the full 
calculation; the dashed curve is obtained by neglecting multiple rescatterings. 



*f°P(P, Q 2 , xjp) PA (b, Zl )p A {b, Z2 )e-™(---)e- (A/2)(1 -^« A 2 d *»^ 



(16) 



The results of the evaluation of the nuclear shadowing correction using Eq. (|16|) at the zero 
impact parameter for anti w-quarks and gluons in 197 Au are presented in Fig. |§] in terms of 
the ratios u A (x, Q 2 , 0)/(AT(0)u n (x, Q 2 , 0)) and g A (x, Q 2 , 0)/(AT(0)g N (x, Q 2 , 0)). The solid 
curves correspond to Q 2 = 4 GeV 2 ; the dashed curves correspond to Q 2 = 10 GeV 2 ; the dot- 
dashed curves correspond to Q 2 = 100 GeV 2 . Note that the factor T(0) = f dzp A {b = 0,z) 
provides the correct normalization of the impulse approximation term, see jsj for details. 
The impact parameter-dependent nPDFs have been tabulated and are available upon request 
from V. Guzey or from the following web cite, http://www.tp2.rub.de/~vadimg/index.html 



They were already used by R. Vogt in the analysis of the J/ip production at RHIC and were 
found to be in a reasonable agreement with the data [43]. 
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FIG. 9: Nuclear shadowing at zero impact parameter: The ratios ua(x, Q 2 , 0)/(AT(0)un(x, Q 2 , 0)) 
and g A (x,Q 2 ,0)/(AT(0)g N (x,Q 2 ,0)) for 197 Au at Q 2 = 4 GeV 2 (solid), Q 2 = 10 GeV 2 (dashed) 
and Q 2 = 100 GeV 2 (dot-dashed). 

V. COMPARISON TO THE DATA AND EVIDENCE FOR HIGHER TWIST EF- 
FECTS 



Our predictions 
to the NMC data 



for the F 2 /(AF 2 ) ratio, where F 2 = (F 2 + F^)/2, can be compared 



44, 



45| . However, since the low-x data points correspond to low Q 2 



we cannot make a direct comparison with those points. Therefore, we simply evaluate 
F 2 A /(AFf) at Q 2 = Q 2 = 4 GeV 2 for the data points with Q 2 < 4 GeV 2 . Figures M and HU 
compare predictions of our leading twist model (upper set of solid curves with the associated 
error bands denoted by dashed curves) to the NMC data on 12 C and 40 Ca |44j; Fig.ll2lmakes 
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FIG. 10: Comparison of the leading twist theory results (upper set of solid curves and associated 
dashed error bands) to the NMC data on F^ a /F^ [44]. The lower set of the solid curves is obtained 
by adding the VMD contribution using Eq. ((T%|) . 

a comparison to the NMC F[ b /(F^) ratio Q] 

One can see from Figs. ^2 and that the agreement between the data points and 
our calculations at low x is poor. Regardless the fact that our model predicts a significant 
nuclear shadowing effect for low-x, x < 10 -3 , see Fig. nuclear shadowing rather rapidly 
decreases when x approaches the values of x probed in fixed-target nuclear DIS experiments. 
Of course, one might argue that we are comparing our predictions at Q 2 = 4 GeV 2 to the 
data with much lower Q 2 values. For instance, in Fig. ^2 for the first five data points, the 
average values of Q 2 are (Q 2 ) = (0.60,0.94,1.4,1.9,2.5) GeV 2 . We have explicitly checked 
that the backward QCD evolution of our predictions down to Q 2 = 2 GeV 2 changes the pre- 
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FIG. 11: Comparison of the leading twist theory results (upper set of solid curves and associated 
dashed error bands) to the NMC data on F^ a /F^ |4J. The lower set of the solid curves is obtained 
by adding the VMD contribution using Eq. ((T%|) . 

dictions only a little. Therefore, since our approach to nuclear shadowing includes the entire 
leading twist contribution to the nuclear shadowing correction (one should keep in mind a 
significant uncertainty due to the unaccounted Reggeon contribution, see Appendix |B]) , the 
disagreement with the NMC low-x data compels us to conclude that the low-x NMC data 
1^1/ could contain significant higher twist effects, which contribute approximately 50% to 
the nuclear shadowing correction to F^. 

At the same time, the A-dependence is reproduced reasonably well (see below) indicating 
that the inadequate modeling of the diffraction at low Q 2 and x is to blame. Indeed, it is very 
natural to have rather significant higher twist effects at small Q 2 since, for this kinematics, 
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FIG. 12: Comparison of the leading twist theory results (upper set of solid curves and associated 
dashed error bands) to the NMC data on Ff b /Ff [45|. The lower set of the solid curves is obtained 
by adding the VMD contribution using Eq. ((T%|) . 

the contribution of small diffractive masses Mx becomes important. Production of small 
diffractive masses Mx is dominated by the production of vector mesons, which is definitely 
a higher twist phenomenon. In this kinematics, the leading twist HI parameterization of 



diffraction 



29| underestimates by approximately factor two the diffractive cross section as 



illustrated by the following estimate. 

Using the definition of the diffractive differential cross section in terms of the diffractive 



structure function F, 



D(3) 



2811 and the relation between the differential cross section on the 



lepton level to the total cross section on the virtual photon level (the Hund convention for 
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the virtual photon flux), cr(7*p — > Xp), the latter can be written as 

atfp - Xp) = fj'° d XlP F 2 Di3 \(3, Q 2 , xjp) . (17) 

Then, if we restrict the integration in Eq. (J 17)) by low diffractive masses Mx, Mj < 1 GeV, 
the resulting cr(7*p — > Xp) can be compared to the cross section of electroproduction of 
vector mesons (dominated by the p meson). For instance, a comparison to the HERMES data 
on exclusive leptoproduction of p° mesons from hydrogen j<46] at low Q 2 and W, (Q 2 ) = 0.83 
GeV 2 and (W) = 5.4 GeV, demonstrates that the calculation using Eq. (|17|) (with the 
restriction Mx < 1 GeV) gives only 40% of the experimental value cr(7*p — > p°p) = 2.04 ± 
0.10 ± 0.43 pb. This observation means that, in the considered case, there is no duality 
between the continuum and resonance contributions to low-mass inclusive diffraction. 

In order to quantitatively study our conclusion about the significant role of the higher 
twist contribution to nuclear shadowing, we explicitly add the contribution of p, <p and to 
vector mesons to our leading twist predictions in the spirit of the vector meson dominance 
(VMD) model. However, since about 50% of the vector meson contribution is already 
contained in the parameterization of inclusive diffraction, we weigh the contribution of the 
vector mesons to nuclear shadowing by the factor 1/2 (this enables us not to double-count the 
vector meson contribution). Therefore, the VMD contribution to the shadowing correction 
to the structure function F 2 reads [l^, \\\ 

SF 2 v j ID 1A-1Q 2 (1-x) ^ a 2 v (Q 2 ) f m 2 ^ 



E 



"V 



AF 2 N 2 2 nFf ^ f 2 \Q 2 + m 2 

P POO POO 

x d 2 b dz x / ^ 2 p A (6, 2l )PA(&^2)cos(A y (^-^))e- (A/2W ^ 2 ^ (M , (18) 
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where Ay = xmy(l + rriy/Q 2 ). The VMD parameters cry and gy assume their usual 
values [l3|. For the nucleon structure function F 2 for low Q 2 and x, we used the NMC 



parameterization |47( for x > 0.006 and the ALLM fit |48( for x < 0.006. In addition, 
since the values of Q 2 , where expression (|THj) is used can be as large as 10 GeV 2 , we take 
into account the effect of the vector meson size decrease by introducing the explicit Q 2 
dependence of cry: cry decreases by the factor of 4 when Q 2 increases from 0.6 — 0.7 GeV 2 
to 10 GeV 2 Q- 

The inclusion of the VMD contribution to the shadowing correction dramatically improves 
the agreement between our calculations and the NMC data for the considered cases of 12 C 
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and 40 Ca, see the lower set of the solid curves in Figs.^Jand^J This supports our conclusion 
about the 50% contribution of higher twist effects to nuclear shadowing in the fixed-target 
nuclear DIS kinematics. 

As seen from Fig.^J the inclusion of the VMD contribution does not significantly improve 
the agreement between our calculations and the data. Since the data points have rather large 
Q 2 , Q 2 > 3.4 GeV 2 , it is natural that the contribution of the p, u and <fi mesons is rather 
small. 

At the same time, the amount of shadowing is sensitive to the diffractive masses up to Q ~ 
1.7 GeV, where an enhancement is also possible as compared to the leading twist fit that we 
employ. One has to emphasize that one is dealing here with xjp > 0.01, where the Pomeron 
component of the diffractive PDFs is known to underestimate the data by a large factor, 
for the recent discussion see 37[. This contribution is usually referred to as the Reggeon 
contribution to diffractive PDFs. However, uncertainties in the theoretical treatment of this 
contribution are very large due to possible interference between the Pomeron and Reggeon 
contribution, unreliable parameter rj for the Reggeon contribution, difficulty to reliably 
extract the Reggeon contribution from the diffractive data. Hence we were forced to neglect 
this contribution in our numerical analysis. Therefore, the theoretical uncertainty due to 
the Reggeon contribution for x > 0.01 is almost out of control and, hence, our calculation 
for x > 0.01 become much less reliable. We refer the reader to Appendix [B] where we 
examine the Reggeon contribution and its influence of our predictions of nuclear shadowing 
and on the comparison of our results to the NMC nuclear DIS data. The domain, where the 
presented leading twist model is best justified, is x < 5 x 10~ 3 and Q 2 > 4 GeV 2 . 

It is worth emphasizing here that our conclusion about the importance of the higher twist 
(HT) effects relies on the existing parameterizations of the leading twist (LT) diffractive 
PDFs which exclude from the fits the region of small diffractive masses. Hence the large /3 
region is, to the large extent, an extrapolation of lower (3 data or effectively due to backward 
evolution of large j3 and high Q 2 data, which is rather not stable. So one cannot a priory 
exclude that there exists a LT parameterization which satisfies a local duality requirement 
in the large (3, low Q 2 region. An investigation of such a possibility would require studying 
diffraction at lower energies than presently accessible at HERA. 

One can also study the A-dependence of nuclear shadowing. Figure presents a com- 
parison of the NMC data on F 2 A /F 2 C Q Q to our model calculations. All data points 
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correspond to x — 0.0125; the F£ ) / F£ data point has Q 2 = 2.3 GeV 2 and all other data 
points have Q 2 = 3.4 GeV 2 . The solid curve is our main prediction, which is obtained as 
a sum of the leading twist shadowing and the VMD contribution evaluated using Eq. (|18|). 
The dashed curve is obtained by increasing the nuclear shadowing correction of the solid 
curve by the factor two. 

As seen from Fig. El our calculation reproduces fairly well the A-dependence of nuclear 
shadowing. However, as was discussed previously in connection with Fig. even with 
the higher twist VMD contribution included, our calculation systematically underestimates 
the absolute value of the shadowing effect. For the heavier nuclei, we should have had a 
2-3 times larger nuclear shadowing, as indicated by the dashed curve. The reason for the 
discrepancy between our calculations and the data is a significant Reggeon contribution to 
hard diffraction at x — 0.0125, which we neglect in our calculation. Because of this, the 
theoretical uncertainty of our results becomes very large for x > 0.01. 

The leading twist model of nuclear shadowing works best in the kinematics where the 
VMD and Reggeon contributions to hard diffraction are small corrections and where the 
data on hard diffraction at HERA were taken, i.e. for Q 2 > 4.5 GeV 2 and x < 5 x 10~ 3 . 
An example of the proper application of the leading twist model of nuclear shadowing is 
presented in Fig. which depicts the A-dependence of / (AF^) at Q 2 = 4 GeV 2 for 
x = 10~ 3 (solid curve), x = 10~ 4 (dashed curve) and x = 10~ 5 (dot-dashed curve). 

Our conclusion about the importance of the higher twist effects at small x and small 
Q 2 in the fixed-target data is in a broad agreement with phenomenological approaches to 
nuclear shadowing, which include both the scaling (leading twist) and lowest mass (p, uj 
and 0) vector meson (higher twist) contributions. (The only possible way to avoid such a 
conclusion would be the local duality scenario which we discussed earlier in this section.) In 



12|, the scaling contribution arises as the effect of the diffractive scattering, quite similarly 



in the spirit to the present work. However, it is difficult to assess the comparability of 
the pre-HERA parameterization of diffraction used in |l2| with the modern HERA data on 
hard diffraction. More importantly, the effect of nuclear shadowing was discussed for the 
nuclear structure function F 2 and not for nPDFs (this comment also applies to all other 
work mentioned below). In other approaches, the scaling contribution results from the qq 
continuum of the virtual photon wave function [13J, or from the contributions of higher 
mass vector mesons or from the aligned qq jets [ill ^|, or from the asymmetric 
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FIG. 13: A-dependence of nuclear shadowing. The NMC data on F^/F^ j3 5 3 are compared to 
our LT+VDM predictions (solid curve). The dashed curve is obtained by scaling up the shadowing 
correction of the solid curve by the factor two. 



qq fluctuations of the virtual photon 18] . Note also that in the case of the real photon 
interaction with nuclei (the accurate data on the real photon diffraction for the relevant 
energies is available, see 5fJ), the shadowing data agree well with the Gribov's theory, see 
the discussion in Ref. [14]. 

We would like to point out that a fairly good description of the NMC data was achieved 
in (^J, which uses the approach to nuclear shadowing based on its relation to diffraction on 
the nucleon. As an input for the their calculation, the authors used the phenomenological 
parameterizations of the inclusive and diffractive structure functions of the nucleon, which 
fit well the inclusive and diffractive data. However, in contrast to our strictly leading twist 
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FIG. 14: The ^4-dependence of nuclear shadowing at Q 2 = 4 GeV 2 . The solid curve corresponds 
to x = 10~ 3 ; the dashed curve corresponds to x = 1CP 4 ; the dot-dashed curve corresponds to 
x = 10" 5 . 



analysis, the phenomenological parameterizations of [2l| have the Q 2 -dependence of the 
form {Q 2 /(Q 2 + where a and b are numerical parameters. Hence, the analysis of 21 1 
effectively includes higher twist contributions, which indirectly confirms our conclusion that 
a good description of the NMC data Q, Q is impossible to achieve without the inc usion 
of the higher twist effects (contribution of vector mesons). 

In addition, we would like to point to other important difference between the present 
analysis and the analysis of [2l( . In order to evaluate nuclear shadowing as a function of Q 2 , 
the authors of apply an equation similar in the spirit to our Eq. (|T4"|) at all Q 2 . As we 
explain in the end of Sect.HIl the application of Eq. (fl4"|) at large Q 2 violates QCD evolution 
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because one then ignores the proper increase of the fluctuations of result of the 

QCD evolution. Also, neglecting proper QCD evolution, one neglects the contribution of 
larger x effects - antishadowing and EMC effects - to the small-x region. The second major 
difference is that the use of the QCD factorization theorem for hard diffraction allowed us to 



make predictions for nPDFs. Since this factorization theorem is not used in [2jj, only nuclear 
structure function F 2 A is considered. Finally, the well-understood and important effect of 
the decrease of the coherence length with the increase of x, i.e. the factor e lXlpmN ^ Zl ~ Z2 \ was 
ignored in 21 1. 

Next we discuss the importance of the next-to-leading order (NLO) effects in the nuclear 
structure function F 2 . Using the LO parameterization for the ratios of the nuclear to proton 
PDFs of Eskola et al. |l| at the initial scale Qq = 1.5 GeV, we perform QCD evolution to 
Q 2 = 10 GeV 2 both with NLO and LO accuracy. The resulting Ff a /(AFg) ratios after 
the NLO evolution (dashed curve) and LO evolution (dot-dashed curve) are presented in 
Fig. E3 For the proton PDFs, we use CTEQ5 parameterizations [42]: CTEQ5M for the 
NLO calculations and CTEQ5L for the LO calculations. For comparison, we also present 
F^ 11 / (AF2) (solid curves and the associated error bands) calculated using our leading twist 
model. The two solid curves correspond to the two scenarios of nuclear shadowing for gluons 
(two models for B g ). A significant difference between the solid and dashed curves in Fig. ITB1 
demonstrates that the effects associated with the NLO QCD evolution and NLO expression 
for the structure function F 2 are important both in the very low-x region and in the x-region 
of the fixed-target data, x > 0.003. This gives us another indication that the LO fits to the 
fixed-target data of must have significant intrinsic uncertainties, especially at low-x, 

where there is no data and the fits are extrapolations. In addition, Fig. demonstrates 
that it is not self-consistent to use the LO fits for nPDFs in the NLO calculations of various 
hard processes with nuclei, which require NLO nPDFs as an input. 



VI. EXTRACTION OF THE NEUTRON F 2n 

The determination of the small- 2 behavior of the valence quark distributions relies heavily 
on the use of the F 2 D and F 2p data for the determination of the F 2p — F 2n difference. The 
main problem is that at x < 0.03, the difference is comparable or even smaller than the 
nuclear shadowing correction for the deuteron, which we denote as 5(x, Q 2 ), 5(x,Q 2 ) = 
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FIG. 15: LO vs. NLO evolution. The LO fits for nuclear PDFs of Eskola ei a/. LU are evolved to 
Q 2 = 10 GeV 2 to LO (dot-dashed curve) and to NLO (dashed curve) accuracy. Our leading twist 
nuclear shadowing predictions are given by the solid curves and the associated error bands. 

1 - F 2 D (x, Q 2 )/(F 2p (x, Q 2 ) + F 2n {x, Q 2 )). Then the F 2p - F 2n difference reads 

F 2p (x, Q 2 ) - F 2n (x, Q 2 ) = 2F 2p (x, Q 2 ) - F 2 (x, Q 2 ) (l + 8(x, Q 2 )) . (19) 

In view of the discussed uncertainties in the model of diffraction (higher twist effects, 
subleading contributions) in the kinematics, where the NMC took their most accurate data 
for fi — D and /i — p scattering, one can hardly use the Gribov theory to calculate reliably 
5(x,Q 2 ). 

However, the same mechanism works both for the deuteron and for heavier nuclei as 
strongly suggested by Fig. EH Hence we can employ the information on the F 2 A /F 2 D ratios 
in order to determine the value of S(x,Q 2 ). Combining the results of the present analysis 
with the results of our analysis of nuclear shadowing corrections to the deuteron F 2 (51 1, 
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we find that 

1/ F 2 c (x,Q 2 )\ F 2 D (x,Q 2 



7 V F 2 D (x,Q 2 )J F 2p (x,Q 2 ) + F 2n (x,Q 2 )- 1 ; 

Using this equation and the NMC nuclear data, we find that 5(0.0125, 2.3 GeV 2 ) ~ 0.011 ± 
0.001. Using the CTEQ5L fit to the nucleon PDFs 42j, we observe that for this kinematics 
nuclear shadowing significantly changes the F 2p (x, Q 2 ) — F 2n (x,Q 2 ) difference 

F 2 Jx, Q 2 ) - F 2n (x, Q 2 ) A + R 

2pK ^ J J^-— i-L = -8^— = -0.63 ± 0.06 , 21 

F 2p (x,Q 2 )-F 2n (x,Q 2 ) l-R 

where R = F 2n / F 2p \ F 2n is the neutron structure function extracted from the deuteron F 2 
ignoring the shadowing correction. In this estimate, we included only errors due to the 
experimental uncertainty in the F 2 / F 2 ratio. 



VII. CONCLUSION AND DISCUSSIONS 



The main results of this paper can be summarized as follows: 

• We explain the derivation of the leading twist theory of nuclear shadowing in DIS on 
nuclei that relates nuclear shadowing in DIS on nuclei to DIS diffraction on the proton. 
The theory enables us to predict nuclear shadowing for individual nuclear PDFs in a 
model-independent way at small x, 10" 5 < x < 10~ 2 . At lar ger x, other nuclear effects 
(antishadowing, EMC effect) and details of the mechanism of diffraction at high xjp 
introduce a large model dependence and uncertainty. 

• Nuclear shadowing corrections to nPDFs are found to be large. In particular, we 
predict larger shadowing than given by the fits by Eskola et al. [l| for gluons for all x 
and for quarks for x < 5 x 10~ 4 . In a stark disagreement with all other approaches, 
we predict larger nuclear shadowing for gluons than for quarks. 

• The presented formalism is applied to evaluate nuclear shadowing for nPDFs at all 
impact parameters. As one decreases the impact parameter, the effect of nuclear 
shadowing increases. 

• The results of our purely leading twist calculations for the F 2 / 'F 2 , Ff°jF 2 N and 
F 2 Ph I F 2 ratios disagree with the corresponding fixed-target NMC data |44|, |45j at low 
x and low Q 2 . While we cannot compare our prediction directly to the data at the 
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Q 2 values of the first five data points, we notice that the backwards QCD evolution 
is small and it does not seem to increase nuclear shadowing. Hence, we conclude that 
the NMC data with Q 2 < 4 GeV 2 are likely to contain a significant amount (about 
50%) of higher twist effects. This is supported by the explicit inclusion of the p, u and 
4> meson contributions to the shadowing correction in the spirit of the vector meson 
dominance model. An alternative scenario would be the existence of a local duality 
pattern for diffraction at xjp > 0.01 where so far no data were taken, see Appendix iBl 
This implies that a leading twist QCD analysis of the low-x and low-Q 2 fixed-target 
data will not produce reliable results for the low-x nuclear PDFs. 

• Using general features of the Gribov theory and the data on A > 2 nuclei, it is 
possible to develop a reliable procedure for the extraction of the neutron F 2n in the 
NMC small-x kinematics. 

• Our predictions for nPDFs, impact parameter-dependent nPDFs and the structure 
function F% for the nuclei of 12 C, 40 Ca, 110 Pd, 197 Au and 206 Pb and for the kine- 
matic range 10~ 5 < x < 1 and 4 < Q 2 < 10, 000 GeV 2 have been tabulated. They 
are available in the form of a simple Fortran program from V. Guzey upon request, 
vadim.guzey@tp2.rub.de. The QCD evolution was carried out using the QCDNUM 
evolution package J3]. 

We must mention that there is a renewed interest in nuclear shadowing because of the 
recent surprising measurements of the suppression of production of hadrons with high p t 
in deuteron-gold collisions at RH1C. It was claimed that the observed suppression is a 
spectacular confirmation of the Colored Glass Condensate model. However, in the RH1C 
kinematics nPDFs are probed at relatively large values of Bjorken x, on average x > 0.01, 
which is beyond the domain of the Colored Glass Condensate model. Since leading twist 
nuclear shadowing is rather weak for x > 0.01, nuclear shadowing cannot be responsible for 
the dramatic effect of the suppression of the hadron spectra at forward rapidities at RH1C, 
see j53l ] - However, shadowing in the forward RHIC kinematics can be observed by selecting 
appropriate two jet production kinematics. 

After the first version of this paper was released, there appeared an analysis j^] which cal- 
culates higher twist effects in shadowing. Similarly to us, the authors come to the conclusion 
that the higher twist effects in the fixed-target kinematics are large. Within uncertainties of 
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Nucleus 


Po (fm -3 ) 


c (fm) 


40 Ca 


0.0039769 


3.6663 


no pd 


0.0014458 


5.308 


197 Au 


0.000808 


6.516 


206p] 3 


0.0007720 


6.6178 



TABLE I: The parameters entering the nuclear one-body density pA( r )- 

their analysis, the higher twist effects could be even responsible for all shadowing observed 
at fixed target energies. So far the connection of the approach of [54] to the Gribov the- 
ory is not clear. In particular, the diagrams, which correspond to vector meson production 
(which dominates the higher twist small x contribution in the Gribov theory, see discussion 
in Sect. EJ) seem to be neglected as a very high twist effect. 
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APPENDIX A: NUCLEAR DENSITY p A 

The nuclear density p A , which enters the calculation of nuclear shadowing in Eq. (|14jh 
was parametrized in a two-parameter Fermi form for 40 Ca, 110 Pd, 197 Au and 206 Pb |5f 



PA{T) = l + exptfr-cVa] ' (A1) 

where r = \J \ b \ 2 + z 2 and a = 0.545 fm and the parameters po and c are presented in TableHJ 
Also note that pa(E,z) was normalized as 2n J °° d\b\ dz\b\pA(b, z) = 1. 
For 12 C, we used 

PA(r)= Po (l + «(£) 2 )e- r2 / a2 (A2) 
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with po = 0.0132, a = 1.403 and a = 1.635 fm. 



APPENDIX B: SUBLEADING (REGGEON) CONTRIBUTION TO NUCLEAR 
SHADOWING 

The analysis of the 1994 HI data on hard diffraction was carried out with the assumption 
that the diffractive structure function F 2 ^ is described by a sum of the effective Pomeron 
(leading) and Reggeon (subleading) contributions 28] 

F 2 D{3 \ X]P , (3, Q 2 ) = fjp /p (xjp) F 2 F ((3, Q 2 ) + fn/pixp) F 2 R (P, Q 2 ) , (Bl) 

where fjp/ p and fm/ p are the so-called Pomeron and Reggeon fluxes; F 2 and F 2 M are the 
Pomeron and Reggeon structure functions. It is important to note that the both terms 



in Eq. (|B1|) are leading twist contributions. In our analysis we use fit 3 of model B of 28 1 
which assumed no interference between the leading and subleading contributions. The latter 
contribution becomes important only for large values of measured xjp, xjp > 0.01. 

In the HI QCD analysis, F 2 M was assumed to be the pion F 2 56| multiplied by a free 
coefficient Cjr to be determined from the data. Unfortunately, the value of Cm is not given 
in the HI publication. Therefore, using Eq. (|B1|) we performed a \ 2 fit to a set of selected 
HI data points with xp > 0.01 and found that Cjr ~ 17. 

The Reggeon contribution to the nuclear shadowing correction to the F 2 structure func- 
tion has the form similar to our master equation ()14j) 



2 |_ 1 + V 2 m 

xF£(f3, Q 2 ) PA (b, Zl) PA (b, Z2 )e ix Frn N (zi-Z2) e -(A/2Kl~i, M )*% f% dz PA (b,z) 



A(A-l) 



v 1 - 



2 r roo poo i>0.1 

d 2 b / dz\ J dz 2 / dxjp 

J~ oo J z\ Jx 



(B2) 



In this equation, the cr^ rescattering cross section is defined as (compare to Eq. (J5J)) 

q2) - F?(M + te r ixF ™ q2) • (B3) 

where it is worth noting the absence of the Reggeon flux since fm./ p {x]p,t = 0) = 1 [28]. 

The theoretical uncertainty associated with the Reggeon contribution originates from the 
uncertainty in the choice of which is the ratio of the real to imaginary parts of the 
subleading exchange amplitude, 

^ + cos7ra K (0) 

Vm = : TTTx — > B4 

sni7ra:jR(0j 
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FIG. 16: The F^-/(AF^) ratio at Q 2 = 4 GeV 2 for Ca-40. The thick solid curves present our main 
result; the thin solid curves present the uncertainty of our predictions; the dashed curve presents 
the fit of p|; the dot-dashed curves is the result of the calculation including both the Pomeron and 
Reggeon contributions to diffraction. The lower dot-dashed curve corresponds to t/jr = —1; the 
upper one corresponds to t/jr = 1. 



where £ = ±1 is the signature factor. The intercept of the Reggeon trajectory, «#?(()), 
was a fit parameter in the HI analysis. The fit to the HI diffractive data gives q;jr(0) 
0.5±0.11(stet)±0.11(s2/s. 



28| . which leads to rjjR = ±1 depending on the signature factor. 
For the p and u meson exchanges, rjm = 1; for the a and / meson exchanges, = — 1. Since 
no attempt was made to separate the contributions with different signatures to the Reggeon 
contribution in the analysis of (inclusive diffraction is not sensitive to the signature of 
the exchange) the value of r/^ is uncertain, — 1 < tjjr < 1. 
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Because of the large |?/2r|, the uncertainty in the choice of r]^ leads to a very significant 
uncertainty in the resulting shadowing correction. This is illustrated in Fig. 1161 presenting 
Ff/(AF?) for Ca-40 at Q 2 = 4 GeV 2 . In Fig. [El the thick solid curves present the 
predictions of our model neglecting the subleading exchange contribution to hard diffraction; 
the thin solid curves present the uncertainty of the predictions; the dashed curve presents 
the result of [j|; the dot-dashed curves is the result of the calculation including both the 
Pomeron and Reggeon contributions to diffraction. The lower dot-dashed curve corresponds 
to 77 jr = — 1; the upper one corresponds to tjjh = 1. As seen from Fig. the variation of 
tjm between its lower and upper limits leads to a dramatic change in the predicted nuclear 
shadowing. 

The influence of the Reggeon contribution on the comparison of our predictions to the 
NMC fixed-target nuclear DIS data is presented in Figs. [T7| and [TH1 The upper set of solid 
curves and the associated dashed error bands is the result of the calculation using only 
the Pomeron contribution to diffraction. The lower set of the solid curves is obtained by 
adding the VMD contribution using Eq. (|18|1. The dot-dashed curves are the results of the 
calculation taking into account the Pomeron, Reggeon and VMD contributions. The lower 
dot-dashed curve corresponds to r)m — — 1; the upper one corresponds to T)m — 1. 

As seen from Figs. El and EE varying rjm in the —1 < t^jr < 1 range, one obtains a wide 
spectrum of predictions for F^ a /F® and F^/F^, which accommodate the NMC data. This 
indicates the possibility of a duality between the higher twist vector meson contribution and 
the leading twist subleading contribution. 
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Hud lint ions in the strength of the interaction arc smaller (see (he discussion in Kef. ||). The 
QCD evolution equations automatically account for the proper increase of the fluctuations of 
the effective cross section around its average value a J s with an increase of Q 2 . This important 
effect is omitted if one attempts to apply Eq. © at Q 2 > Qq with a Q 2 -dependent a^^Q 2 ). 
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